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ABSTRACT 



The postburst evolution of fireballs that produce 7-ray bursts is studied, assuming the ex- 
pansion of fireballs to be adiabatic and relativistic. Numerical results as well as an approximate 
analytic solution for the evolution are presented. Due to adoption of a new relation among 
t, R and 7 (see the text), our results differ markedly from the previous studies. Synchrotron 
radiation from the shocked interstellar medium is attentively calculated, using a convenient set 
of equations. The observed X-ray flux of GRB afterglows can be reproduced easily. Although 
the optical afterglows seem much more complicated, our results can still present a rather satis- 
factory approach to observations. It is also found that the expansion will no longer be highly 
relativistic about 4 days after the main GRB. We thus suggest that the marginally relativistic 
phase of the expansion should be investigated so as to check the afterglows observed a week or 
more later. 
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1 INTRODUCTION 



Since their discovery nearly thirty years ago (Klebesadel, Strong &; Olson 1973), 7-ray bursts 
(GRBs) have made one of the biggest mysteries in astrophysics (Fishman k. Meegan 1995), 
primarily because they have remained invisible at longer wavelengths. The situation started to 
change drastically in 1997 due to the Italian-Dutch BeppoSAX satellite (Piro et al. 1995), whose 
prominent observations led to the discoveries of multi- wavelength counterparts of several GRBs: 
GRB970111 (Costa et al. 1997a), GRB970228 (Costa et al. 1997b), GRB970402 (Feroci et al. 
1997; HcisG et al. 1997) and GRB970508 (Costa et al. 1997d). The corresponding afterglows in 
X-rays (GRB970228: Costa et al. 1997c,e; GRB970402: Piro et al. 1997a; GRB970508: Piro et 
al. 1997b), in optical band (GRB970228: Groot et al. 1997, van Paradijs et al. 1997, Sahu et al. 
1997, and Galama et al. 1997; GRB970508: Bond 1997), and in radio band (GRB970508: Frail 
et al. 1997) were observed with unprecedented enthusiasm and collaboration, and the results 
are exciting. GRB970228 seems to be associated with a faint galaxy (van Paradijs et al. 1997), 
and the redshift of the optical counterpart of GRB970508 was even determined to be between 
0.835 and 2 (Metzger et al. 1997). Very recently it is also reported that X-ray afterglows 
of GRB970616 (Marshall et al. 1997a; Murakami ct al. 1997a), GRB970828 (Remillard et 
al. 1997) and possibly GRB970815 (lAU Circular No. 6718) were observed due to efficacious 
cooperations between Compton Gamma-Ray Observatory and Rossi X-ray Timing Explorer. 
These observations strongly suggest that GRBs originate from cosmological distances. The 
fireball model (Goodman 1986; Paczyhski 1986; Rees & Meszaros 1992,1994; Meszaros & Rees 
1992; Meszaros, Laguna & Rees 1993; Meszaros, Rees & Papathanassiou 1994; Katz 1994; Sari, 
Narayan & Piran 1996) becomes the most popular and successful model for GRBs, although 
other models, such as a hypernova scenario proposed by Paczyhski (1997) can not be eliminated 
now. 

After producing the main GRB, the cooling fireball is expected to expand as a thin shell 
into the interstellar medium (ISM) and generate a rclativistic blastwave, although whether the 
expansion is highly radiative (Vietri 1997a,b) or adiabatic is still controversial. Afterglows at 
longer wavelengths are produced by the shocked interstellar medium (ISM) (Meszaros k, Rees 
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1997). Much analytical work on GRB afterglows has been done (Waxman 1997a, b; Tavani 
1997; Sari 1997; Meszaros, Rees & Wijers 1997; Dai & Lu 1997), and it has been found that for 
adiabatic expansion, R oc oc t~^/^, where R is the shock radius measured in the burster's 

static frame, 7 is the Lorentz factor of the shocked ISM measured in the observer's frame and 
t is the observed time. These scaling laws are valid only at the ultrarelativistic expansion stage 
(7 » 1)- 

The purpose of this work is to study numerically the evolution of adiabatic fireballs ap- 
propriate from the ultrarelativistic expansion stage to the mildly relativistic expansion stage. 
We show that although the scaling law between 7 and t at the ultrarelativistic expansion stage 
obtained in this work is the same as above, our coefficient for 7 differs dramatically from that of 
Waxman (1997a, b). A detailed set of equations are presented to calculate synchrotron radiation 
from the shocked ISM. We see that radiation during the mildly relativistic phase (2 < 7 < 5) 
which was obviously neglected in the previous studies is quite important. 

The structure of this paper is as follows. In Section 2 we calculate the dynamical evolution 
of the relativistic shock. The difference between our results and those of Waxman is stressed. 
Synchrotron radiation from the shocked ISM is formulated and compared with observations in 
Section 3, and a brief discussion is given in the final section. 

2 DYNAMICAL EVOLUTION OF FIREBALLS 

A fireball with total initial energy Eq and initial bulk Lorentz factor r) = Eq/Mqc^, where Mq 
is the initial baryon mass and c the velocity of light, is expected to radiate half of its energy in 
7-rays during the GRB phase (Sari & Piran 1995), either due to an internal-shock or an external- 
shock mechanism (Paczyriski &; Xu 1994; Vietri 1997b, and references therein). Subsequently 
the fireball will continue to expand as a thin shell into the ISM, generating an ultrarelativistic 
shock. 
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The jump conditions for the shock can be described as (Blandford & McKee 1976): 
n' = (47 + 3)n, 



(1) 



e' = (47 + 2))^nmp(? 



(2) 



(7 + l)(47-lf 
87 + 10 



(3) 



where n is the number density of the unshocked ISM, n' and e' are the number density and 
energy density of the shocked ISM in the frame comoving with the shell, mp is the proton mass, 
and r is the Lorentz factor of the shock in the observer's frame. These equations are valid 
for describing ultrarelativistic shocks as well as mildly relativistic shocks (Blandford & McKee 
1976). For 7 » 1, we have e' = A'y'^nMpC? and T = y/2'j. We will assume that the shocked ISM 
in the shell is homogeneous. 

As usual, the expansion of the fireball is thought to be adiabatic, during which energy is 
conserved, so we have 



where Vg is the observed velocity of the shock and R{1 — Vs/c)j is the thickness of the shocked 
ISM in the comoving frame. / is defined as / = 47ri?'^nmp/(3Mo). Using equation (1), equation 
(4) can be further expressed as: 




(4) 




(5) 



For J ^ 1, equation (5) generates an approximate relation: 87^/ « r//2, or equally. 



7^i?^ Mori/{8Trnmp) = Eo/{8TmmpC^). 



(6) 



In order to study the evolution of 7, we should add two equations. Photons observed at a 
time interval of At are in fact emitted at an interval of Ati,: 



where tb is the time measured in the burster's static frame and v is the velocity of the shocked 
ISM in the observer's frame, while Atj, and variation of radius (Ai?) are related due to 

Vr^ - 1 

AR = ^ cAtb. (8) 

Equations (2), (3), (5), (7) and (8) present a perfect description of propagation of the 
shock. Given initial values of Eq and Mq, R{t) and j{t) can be evaluated numerically. We 
take Eq = 10^^ ergs, n = lcm~^, and Mq = 10~^Mq and 2 x 10~^Mq. Using these equations, 
we numerically study the evolution of 7 and R as functions of observed time. Our results are 
plotted in Figs. 1 and 2, where only the dot-dashed line corresponds to Mq = 10~^Mq, and the 
other lines to Mq = 2 x 10'^ Mq. It can be seen that Mq only influences the early-time evolution 
of the fireball. 

Under the assumption that 7 S> 1, we can derive a simple analytic solution, 

i?^ fii i?^ + Akt, (9) 

(10) 

where k = Eq / {Aimmpc) , and Rq is the initial value of R. In Figs. 1 and 2, we also plot the 
results of equations (9) and (10). It is clear that at early time, when 7 3> 1, the analytic solution 
is a quite good approximation for the numerical results. 

If S> -Ro) that is t S> r, where r refers to the duration of the main GRB, then we can 
rewrite equations (9) and (10) as: 

R « 8.93 X lQ^^El{^n~^'H'^/^cm, (11) 

7 PS l92>El{^n^^'^t-^/^, (12) 

where Eq = IQ^^E^i ergs, n = Ini cm and t is in units of second. We find that when 7 3> 1 
and T, equations (11) and (12) also fit numerical results very well. Equations (11) and (12) 
are scaling laws for R and 7 mentioned in Introduction. 
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We next compare equation (12) with the previous studies. The scahng law for 7 is quite 
simple, but the coefficient should be treated with great carefulness, since it may affect syn- 
chrotron radiation and observational behaviors severely (Sari 1997). Waxman (1997a,b) derived 
a result with a larger coefficient: 7 ^ 332El{^n^ ^^^t^^/^. As shown in the next section, this will 
result in very strong radiation in both X-ray and optical bands. We have plotted his results for 
R{t) and 7(t) in Figs. 1 and 2. The discrepancy is noticeable. The difference between his and 
ours is due to the fact that he has used a relation among t, R and 7: t = R/ (27^C), which may 
be incorrect for the evolution of ultra-relativistic adiabatic fireballs (Sari 1997), and wc have 
adopted another relation: t = i?/(87^c), which is obtained by integrating equations (7) and (8). 

We want to emphasize that the equations in this section are correct only for ultrarelativistic 
blastwaves (7 2> 1) and mildly relativistic blastwaves (2 < 7 < 5). This means that after about 
three days our numerical results for the parameters used in the above might be spurious. Of 
course, our calculation can be considerably prolonged by adjusting the parameter (E^i/ni). 



3 X-RAY AND OPTICAL AFTERGLOWS 

3.1 Synchrotron radiation 

Electrons in the shocked ISM are highly relativistic. Inverse Compton cooling of the electrons 
may not contribute to emission in X-ray and optical bands we are interested in. We will consider 
only synchrotron radiation below. The electron distribution in the shocked ISM is assumed to 
be a power-law function of electron energy, as expected for shock acceleration, 

-J^ ^ le^^ ilmin < 7e < Imax), (13) 

"7e 

where ^ymin and ^ymax are the minimum and maximum Lorentz factors, and p is the index varying 
between 2 and 3. We suppose that the magnetic field energy density (in the comoving frame) is 
a fraction of the energy density, B'^/Stt = ^%e', and that the electrons carry a fraction of 
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the energy. 7^i„ is determined by 



rup \ [p — 2 



We estimate ^max by equating, as usual, the electron acceleration timescale with the synchrotron 
cooling timescale, and find (Meszaros, Laguna & Rees 1993; Vietri 1997a) 

^max ~ 10«S'-V2. (15) 

The spectral property of synchrotron radiation is clear (Rybicki k, Lightman 1979). In the 
comoving frame, the characteristic photon frequency is = e-B'7^j„/(27rmec), where e is the 
electron charge. The spectral peaks at fmax ~ 0.29^^. For frequency u S> Vmax^ the flux density 
scales as 5,^ oc where a = (jp — l)/2, and for v <^ i^max, Si, oc i/^/^. Below we formulate 
synchrotron radiation. First, using equation (1), we further express equation (13) as 

g = C'.- (16) 



where 



G' = b-l)7^7n(47 + 3)n. (17) 



Second, the synchrotron radiation power emitted per unit volume is 

'Ymin 

with 

+ 0O 

F{x) = x I K5/3{t)dt, (19) 



and 
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where i^5/3(i) is the Bessel function. The specific intensity at frequency v in the comoving frame 
is thus written as 

h,co = - ^)7- (21) 

The observed frequency and specific intensity are related to u and ij^^co due to the 

following relativistic translations (Rybicki & Lightman 1979): 

i^e = (1 + v/c)-fu, (22) 
I,^,e = (l + ^/c)VW (23) 

In previous researches, it is customary in analytic derivations to consider that at a given 
time t the emitting surface is located at 2j^ct and that the disk seen by the observer has a radius 
7ct, as it would be in the absence of the deceleration (i.e., an ellipsoid). Then the relativistic 
fireball expanding isotropically will produce an observed flux S,y^^^ = 7r(7ci)^/,y^,®/I?^, where 
D is the distance from the source to the Earth (Rees 1966). However, the deceleration dynamics 
was recently investigated in more detail and it was found that due to the deceleration the 
emitting surfaces become distorted ellipsoids, and at sufficiently late times, most of the light 
(either bolometric or in a given band) comes from a ring-like surface (Waxman 1997c; Panaitescu 
&: Meszaros 1997). For a given observed frequency band, as the shocked fluid is decelerated, 
the peak frequency {vpeak = li^max) crosses the observed band, and the region radiating in that 
band shrinks from the full disk to a narrow ring. According to Panaitescu &; Meszaros (1997), 
at energies far below or above Upeak, the ratios of the equivalent radii of the emitting surfaces to 
7ci are approximately constant in time. For example, these ratios are 2.8 and 3.3 respectively, 
by assuming a homogeneous ISM and an adiabatic expansion. We use the following equation to 
evaluate the observed flux density, 

_ 7r(fc7ct)^J^g„e 
'^'^e>® - ^2 ' 

where A; = 2.8 for i^® < Vpeak A; = 3.3 for > Vpeak- Since Vpeak enters X-ray and optical 
bands very quickly, the "visible" zone acts as a narrow ring for most of time and Equation (24) 



should be accurate enough for our calculations. So we get the observed flux density 



4(1 + ^)^(1-7)^' 



Q R{kct)^j{uM) 

£>2 



(25) 



This equation shows that iSjyg,© strongly depends on 7, so that the coefficient in Equation (12) 
should be treated carefully (Sari 1997). The observed flux by a detector is an integral of Si,^^,^: 



where f„ and are the upper and lower frequency limits of the detector. 
3.2 Comparison with observations 

Following the numerical solution in Section 2, we continue to calculate the afterglow in X-ray 
and optical bands, using equation (26). Some of the parameters are taken as follows: p = 2.5, 
= 0.5, and = O-l, which are required by the spectral and temporal properties of GRBs 
(Sari et al. 1997; Wijers, Rees & Meszaros 1997). The distance D is set to be 3 Gpc. In order to 
get an X-ray flux, we integrate equation (26) from 0.1 keV to 10 keV, since it is approximately 
in the bands available for BeppoSAX and ASCA. As for optical flux, we use Equation (25) to 
calculate the R band flux density {Sr). Our numerical results for X-ray flux (Fx) and Sr are 
illustrated in Fig. 3 and Fig. 4. 



Also plotted in Figs. 3 and 4 are some observed data, which would make it possible for us 
to see to what extent the model could agree with observations. The X-ray data are quoted from: 
(1) Wijers, Rees & Meszaros (1997); (2) Frontera et al. (1997); (3) Costa et al. (1997a); (4) 
Butler et al. (1997); (5) Feroci et al. (1997); (6) Piro et al. (1997a); (7) Costa et al. (1997d); 
(8) Piro et al. (1997b); (9) Marshall et al. (1997); (10) Murakami et al. (1997a); (11) Remillard 
et al. (1997); (12) Marshall et al. (1997b); (13) Murakami et al. (1997b); and (14) Greiner 
et al. (1997). Please note that since different detectors work in different bands and here we 
have converted their flux data into 0.1 — 10 keV band linearly, errors to a factor of two are thus 
possible. In Fig. 4 the observed R band flux densities are quoted from: (1) Wijers, Rees & 
Meszaros (1997); (2) Galama et al. (1997); (3) Fruchter et al. (1997). 




(26) 
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Although the observed GRBs are expected to reside at different distances and their intrinsic 
parameters such as Eg, n, p, and may vary markedly, it can be seen from Fig. 3 that the 
observed X-ray data are really quite easy to be reproduced by our model. The optical afterglow 
from GRB 970228 can be fitted quite well for t < 3 days. However after about three days the 
theoretical light curve shows too sharp a decline. As mentioned in Section 2, this might be 
spurious, since our calculations are reliable only before about three days. Because the X-ray 
afterglow of GRB970228 was observed 11.6 days later and the optical afterglow was observed 
more than 6 months later, we suggest that a marginally relativistic (1.0005 < 7 < 2, or 10^ 
km/s < u < 2.6 X 10^ km/s) model should be considered so as to provide a perfect description 
for the afterglow. We noticed that Wijers et al.(1997) have pointed out that when the GRB 
remnant becomes nonrelativistic and enters the Sedov- Taylor phase, the optical flux will decline 
as Fqp oc i(3-i5a)/5^ decaying faster than earlier time. 

On 4 September the Hubble Space Telescope observed the afterglow of GRB 970228 for the 
third time and found that the optical transient has faded to V ^ 28"*. (Pruchter et al. 1997), 
which corresponds to the R-band magnitude ^ 27.0, being well consistent with the power- 
law extrapolation of earlier data. This suggests that the radiation during the whole period 
may be emitted through one mechanism. Although this seems to have confirmed the fireball 
model, wc argue below that it may be a puzzle. If supposing 7 dacays as j{t) = (200 ~ 
300)(£^5i/ni)^/^t~^/®, the relativistic condition 7 > 2 for i = 6.25 months will require (£^5i/ni) > 
(10^ ~ 10^), which seems quite unlikely. This difficulty may be overcome by assuming that the 
ISM is not homogeneous, so that the shock can keep to be relativistic for more than 6 months. 
Another possibility is that 7 does fall below 2 several days after the main GRB, but the radiation 
during the marginally relativistic phase could account for the long-term optical afterglow. This 
point should be investigated in more details. 

The optical afterglow of GRB 970508 shows a first-rising-then-decreasing behavior (Castro- 
Tirado et al. 1997; Vietri 1997b; Djorgovski et al. 1997). The hght curve peaks at about t = 2 
days. This behavior can not be explained by the simple shock model described here, which only 
predicts a peak at t = 10^ ~ 10^ s. We suggest that an inhomogeneous ISM with some clumps 
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might account for it. 



4 DISCUSSION 

BeppoSAX has made an important breakthrough in GRB researches. GRBs are widely believed 
to be produced by relativistically expanding blastwaves, or namely fireballs, at cosmological 
distances. Catastrophic events such as merging of neutron star binaries (Narayan et al. 1992; 
Vietri 1996), failed supernovae (Woosley 1993), collapse of magnetized white dwarfs (Usov 1992), 
accretion-induced phase transitions of accreting neutron stars (Cheng & Dai 1996), and hyper- 
novae (Paczyhski 1997) have been suggested as possible cosmological models of GRBs. One 
expects that extensive observational and theoretical investigations on GRB afterglows should 
be helpful to providing much more important clues to studies of the GRB origin. 

We have shown in this paper that a relativistic fireball expanding adiabatically into the 
uniform interstellar medium can roughly explain the afterglows of five observed GRBs, especially 
their X-ray flux. We have made a detailed numerical study of the expansion and derived an 
approximate analytic solution. We would like to stress that the relation among t, R, and 7 
should he t = R/{8j'^c) or At = AR/{2^'^c). Our results indicate a smaller coefficient for j{t): 
7(t) 193(£;5i/ni)V8t-3/8^ ^Y^^^Yi differs noticeably from that of Waxman's (1997a,b). The 
difference in the coefficient results in great differences in observational effects. We present a 
set of equations to describe the synchrotron radiation from the shocked ISM, and calculate the 
X-ray and optical flux. It is found that an expanding fireball will no longer be highly relativistic 
(7 < 2) about 3 to 4 days after the main GRB. This leads to our suggestion that a marginally 
relativistic expansion model analogous to (but still much faster than) that for a supernova 
(Woosley & Weaver 1986) should be established. 

We would like to thank Dr. Ralph Wijers for many good suggestions and helpful comments 
that have led to an overall improvement of this paper. This work was supported by the National 
Natural Science Foundation, the National Climbing Programme on Fundamental Researches 
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Figure Captions 



Figure 1. Evolution of the shock radius. The sohd hne is our numerical result and the dashed 
line our approximate analytic solution (equation [9]). Waxman's result is plotted by the dotted 
line. Only the dot-dashed line corresponds to Mq = 10~^ Mq, and the other lines to Mq = 
2 X IQ-^Mq. 

Figure 2. Evolution of 7. The solid line is our numerical result and the dashed line our 
approximate analytic solution (equation [10]). Waxman's result is plotted by the dotted line. 
Only the dot-dashed line corresponds to Mq = W~^Mq, and the other lines to Mq = 2x 10~^Mq. 

Figure 3. 0.1—10 keV flux {Fx) vs. time. Fx is in units of ergs cm~^ s~^. The full line is our 
numerical result. For the observed fluxes, please see the text. 

Figure 4. Optical flux density (Sr) vs. time. Sr is in units of ergs cm" -2 Hz-^ The solid 
line is our numerical result. For the observed flux densities, please see the text. 



17 




(uio)jj §oi 



